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ABSTRACT 

We augment the two widest /deepest 1.4 GHz radio surveys: the NRAO VLA Sky Survey (NVSS) 
and the Faint Images of the Radio Sky at Twenty-Centimeters (FIRST), with the mean epoch in 
which each source was observed. We use these catalogs to search for unresolved sources which vary 
between the FIRST and NVSS epochs. We find 43 variable sources (0.1% of the sources) which vary 
by more than 4cr, and we construct the mean structure function of these objects. This enables us to 
explore radio variability on time scales between several months and about five years. We find that on 
these time scales, the mean structure function of the variable sources is consistent with a flat structure 
function. A plausible explanation to these observations is that a large fraction of the variability at 
1.4 GHz is induced by scintillations in the interstellar medium, rather than by intrinsic variability. 
Finally, for a sub sample of the variables for which the redshift is available, we do not find strong 
evidence for a correlation between the variability amplitude and the source redshift. 
Subject headings: radio continuum: general — ISM: general — quasars: general 



1. INTRODUCTION 

Variability of radio sources at low frequencies is mainly 
attributed to propagation effects (scintillations) induced 
by large scale electron density inhomogeneities in the 
Inter-Stellar Medium (ISM; e.g., Hunstead 1972; Rick- 
ett, Coles & Bourgois 1984; Rickett 1990; Ghosh & Rao 
1992). The predicted variability structure function of 
compact radio sources due to scintillations (e.g., Bland- 
ford & Narayan 1985; Goodman & Narayan 1985; Bland- 
ford et al. 1986; Hjellming & Narayan 1986) is roughly 
consistent with the typically observed structure function, 
at least below 5 GHz (e.g., Qian et al. 1995; Gaensler & 
Hunstead 2000). These models predict a rise in the struc- 
ture function up to time scales of ~ 10 days at rs 5 GHz 
and up to ~ 100 days at » 0.5 GHz, followed by a flat- 
tening of the structure function. 

Specifically, Qian et al. (1995) analyzed radio observa- 
tions of the compact radio source 1741—038 (z = 1.054) 
taken in several frequencies between 1.5 and 22 GHz. 
They compared the observed structure functions with 
theoretical models for scattering by an extended Galac- 
tic medium, with and without a thin screen component. 
They reported that for frequencies below about 5 GHz 
the observations are consistent with a scattering by an 
extended Galactic medium and a thin screen. However, 
above this frequency they found excess variability rela- 
tive to the models. Moreover, at these high frequencies, 
the structure function continues to rise towards longer 
time scales. They suggested that at high frequencies 
(> 5 GHz) some of the variability of this radio source 
is intrinsic. This general picture is also supported by 
Mitchell et al. (1994). 

Gaensler & Hunstead (2000) studied the variabil- 
ity of 55 radio calibrators observed by the Molonglo 
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Observatory Synthesis Telescope (MOST) at 843 MHz. 
They constructed the structure function for 18 variable 
sources. For the majority of these variable objects the 
structure function flattens on time scales of a few hun- 
dreds days. Furthermore, they confirmed early results 
(Condon et al. 1979) which found that the variability 
amplitude is increasing as a function of the source spec- 
tral index a, defined by f v oc v a , where f v is the specific 
flux at frequency v. This is attributed to a correlation 
between the source angular size and spectral index. An- 
other confirmation for the importance of Galactic scin- 
tillations is that radio variability depends on Galactic 
latitude (e.g., Spangler et al. 1989; Ghosh & Rao 1992; 
Gaensler & Hunstead 2000). We note however that Rys 
& Machalski (1990) did not find evidence for increasing 
fraction of 1.4 GHz variability for sources brighter than 
100 mJy at low Galactic latitudes. 

Lovell et al. (2008) presented results from the 
Micro- Arcsecond Scintillation-Induced Variability (MA- 
SIV) survey conducted at 5 GHz. Among their findings: 
half of the sources they monitored exhibit 2%-10% rms 
variations on time scales over two days. They also found 
that the structure function of the variable sources rises 
on time scales of a few days, and the variability ampli- 
tudes correlates with the Ha emission at the direction 
of the sources. Furthermore, there is evidence that the 
variability amplitude decrease with redshift above z » 2, 
presumably due to evolution of the source size with red- 
shift (see however Lazio et al. 2008). 

Here we compare the two widest/deepest 1.4 GHz sky 
surveys taken using the Very Large ArrajQ (VLA), and 
search for variable sources. We use these datasets to 
construct the average structure function on time scales 
between several months and about five years and to look 
for indications for intrinsic variability of these sources. 

4 The Very Large Array is operated by the National Radio As- 
tronomy Observatory (NRAO), a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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In we present augmented versions of the FIRST and 
NVSS catalogs which contains the mean time in which 
each source was observed. In $3] we cross correlate the 
two catalogs, and in §4] we construct the structure func- 
tion of the variable sources. Finally, we discuss the re- 
sults in [JS] 

2. THE CATALOGS 

The NVSS observations were carried out between June 
1993 and April 1999, while the FIRST survey observa- 
tions were conducted between March 1993 and Septem- 
ber 2002. Therefore, these observations provide a long 
baseline to the structure function analysis. Constructing 
the structure function of variable sources requires knowl- 
edge of their fluxes at multiple epochs and the time at 
which the observations were taken. However, the NRAO 
VLA Sky Survey (NVSS; Condon et al. 1998) and the 
Faint Images of the Radio Sky at Twenty-Centimeters 
(FIRST; Becker et al. 1995) source catalogs do not con- 
tain the time at which each source was observed. The 
reason for this is that the observing times are not well 
defined. Images in both surveys were taken by scanning 
the sky in a hexagonal grid in which observing points are 
separated by 26'. Both surveys were obtained using the 
VLA, in which the full width at half power at 1.4 GHz is 
= 31'. Each primary beam field of view was truncated to 
radii of 24' and 23.5' for the NVSS and FIRST surveys, 
respectively. Therefore, each point on the sky effectively 
contains information from roughly four different point- 
ings taken at different times. 

In order to obtain the epoch at which each source was 
observed we downloaded from the VLA archiv^ the list 
of all observing scans which are associated with each sky 
survej0. 

Next, we cross-correlate the list of observing scans for 
each project with its catalog of sources. We use a match- 
ing radii equal to the truncation radii of 24' and 23.5' 
for the NVSS and FIRSlQ surveys, respectively. This 
enable us to estimate for each source, the number of ob- 
servations (iV bs), its mean observing time (over all scan 
mid times, t cp ) and the time span within which these 
observation were obtained (5t) . 

The products are versions of the FIRST and NVSS cat- 
alogs that contains the observing time, number of obser- 
vations, and time span of observations for each sourc<fl 
In Tables HH2] we present a version of these catalogs con- 
taining the source coordinates and the observing time 
information for each source. 

Figure [1] presents histograms of the number of individ- 
ual snapshots used to compose each source image. For 
both surveys, the typical number of snapshots per com- 
posite image is 3-4 (see also Helfand et al. 1996). In 
Figure [5] we show histograms of St for the two surveys. 
This figure suggests that most of the images are made 
from snapshots taken within a few weeks range of each 
other. 



TABLE 1 

Observing times of NVSS sources 



J2000 RA 

deg 


J2000 Dec 

deg 


/ 
mJy 


Of 

mJy 


N obB 


Eep 
days 


St 
days 


194.89704 


-40.37908 


3.60 


0.70 


1 


0220.698 


0.000 


249.00512 


-40.37366 


99.70 


3.60 


1 


0220.840 


0.000 


199.42008 


-40.36441 


24.30 


1.20 


1 


0220.699 


0.000 


248.90470 


-40.36169 


3.00 


0.60 


1 


0220.840 


0.000 


212.86787 


-40.35236 


2.70 


0.50 


1 


0220.737 


0.000 



Note. — A version of the NVSS catalog containing the source posi- 
tion, mean observing time (t ep ), number of scans (iV bs), and the time 
span over which the scans were obtained (St). The observing time, £ ep , 
is given in JD — 2450000 days, where JD is the Julian day. / is the peak 
flux density and af is the error in the peak flux density. We note that 
our cross correlation of the observing scans and the sources catalog was 
not able to produce the observing times of 591 sources. This table is 
published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion of the full table is shown here for guidance regard- 
ing its form and content. The table is sorted by Declination, therefore 
the sources listed here are near the edge of the survey footprint and 
have a single observation. 



TABLE 2 

Observing times of FIRST sources 



J2000 RA 


J2000 Dec 


/ 


Of 


N obs 


tep 


St 


deg 


deg 


mJy 


mJy 




days 


days 


354.74904 


-11.39882 


1.68 


0.14 


1 


0593.087 


0.000 


5.19669 


-11.39761 


1.02 


0.15 


1 


0593.118 


0.000 


6.82096 


-11.39611 


1.11 


0.15 


1 


0595.073 


0.000 


359.90417 


-11.39585 


1.39 


0.15 


1 


0593.102 


0.000 


6.70133 


-11.39458 


1.17 


0.14 


1 


0595.073 


0.000 



Note. — Like Table[T]but for the FIRST catalog. We note that our 
cross correlation of the observing scans and the sources catalog was not 
able to produce the observing times of 379 sources. 
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Fig. 1. — Histogram of number of observations from which each 
source image was composed. Bin size is 1 observation. The upper 
panel is for NVSS sources, while the lower panel is for FIRST 
sources. 



5 https:/ /archive. nrao.edu/archivc/advqucry.jsp 

6 These are observing code AC308 for the NVSS catalog and 
AB628, AB879 and AB950 for the FIRST catalog. 

7 We use FIRST catalog version 20080716. 

8 These are approximate observing times since we do not know 
if all the data was used in the reduction process of the FIRST and 
NVSS. 



We note that a transient search based on the compar- 
ison of the FIRST and NVSS catalogs was presented in 
Levinson et al. (2002) and discussed in Gal- Yam et al. 
(2006) and Ofek et al. (2010). However, these previous 
efforts did not use the observing time of the sources. 
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Fig. 2. — Histogram of the time span of observations from which 
each source image was composed (St). Bin size is 1 day. Panels like 
in Fig.[l] The last bin on the right-hand side is for all St > 59 days. 



3. CROSS CORRELATION OF THE FIRST AND 
NVSS CATALOGS 

The FIRST catalog contains 816,331 sources, brighter 
than about 1 mJy, mainly in the North Galactic cap. 
About 81% of these sources have St < 30 days. The 
NVSS catalog contains 1,773,484 objects with 6 > 
— 30deg, brighter than about 3.5 mJy, of which = 67% 
have St < 30 days. We select all the FIRST sources which 
deconvolved major and minor axes equal zero (i.e., point 
sources), the peak flux density is larger than 5mJy, all 
the scans composing their flux measurement (St) were 
taken within 30 days, and which are isolated from any 
other FIRST source (of any kind) within 60". Since the 
resolution of the NVSS is 45", nine times coarser than 
that of the FIRST survey, the last step is designed to re- 
move NVSS sources which flux may be contaminated by 
multiple FIRST objects. We found 6463 FIRST sources 
that satisfy these criteria. Next, from the NVSS catalog 
we select all sources with St < 30 days - we find 1,183,620 
sources that satisfy this criterion. 

Then, for each object in the subset of the FIRST cata- 
log we search for a source in the subset of the NVSS cata- 
log which is found within 15" of the FIRST objec10. We 
found 4367 matched sources. These matches represent 
point sources for which we have both a FIRST and NVSS 
flux-density measurements. We note that only 68% of the 
FIRST sources in this list have NVSS matches. This is 
mostly due to the fact that we used only NVSS sources 
with St < 30 days. 

Next, we would like to compare the fluxes of unresolved 
NVSS and FIRST sources. However, systematic biases 
in the NVSS and/or FIRST flux calibration could effect 
our analysis. Condon et al. (1998) and Becker et al. 
(1995) discussed photometric errors such as the CLEAN 
bias, and they made corrections to their source catalogs 
flux densities. For both the NVSS and FIRST catalogs 
the CLEAN bias is « —0.3 mJy, i.e. of order the rms 
noise in the images. There is also a well-known discrep- 

9 The median astrometric error for 5 mJy sources in the NVSS 
catalog is about 3", and 99.4% of the errors of such sources are 
smaller than 15" 



1.005 



0.995 



0.99 



0.985 



0.98 



0.975 



10 12 14 16 

Flux density [mJy] 



20 



Fig. 3. — The mean of the flux ratio between matched individual 
NVSS and FIRST point sources, as a function of flux density (black 
circles) . The horizontal "error bars" represent the bin size in which 
the mean flux ratio was measured. The expected amplitude of the 
Eddington-like bias (see text) is shown in gray line. 

ancy between integrated fluxes of extended sources in 
FIRST and NVSS, owing to resolution effects (Blake k 
Wall 2002). In Figure |3] we show the mean of the peak- 
flux ratio between matched individual NVSS and FIRST 
point sources, as a function of flux density (black circles). 

This figure suggests that at the faint end, NVSS fluxes 
are systematically lower than FIRST fluxes. Blake & 
Wall (2002) already reported this effect, although with 
an opposite direction and higher amplitude. However, 
Blake & Wall (2002) looked at both resolved and unre- 
solved sourcetQ while we are interested only in point 
sources. Probably the most important reason for this 
trend is related to the fact that the NVSS and FIRST 
surveys have different resolutions. Another secondary ef- 
fect is a bias similar to the Eddington bias (Eddington 
1913). This is because sources which flux expectancy 
value is smaller than 5 mJy (our FIRST flux cut) and are 
detected in the FIRST survey above 5 mJy (due to mea- 
surement errors) have > 50% chance to have an NVSS 
flux density below 5 mJy. This effect is amplified by the 
fact that faint sources are more common, per unit flux 
density, than bright sources. The estimated amplitude 
of this bias, based on simulations, is shown by the gray 
line in Figure [3J 

Given these results we correct the NVSS fluxes of 
the matched point source^]] by the amount interpolated 
from the black circles in Figure [3] Above 20 mJy we 
assume that the correction factor is 1 . We neglect the ef- 
fect of the Eddington bias, since its expected amplitude 
is negligible (see Fig. [3]). We note that the maximum 
amplitude of the applied correction is smaller than the 
rms noise in the NVSS measurements and comparable 
to rms noise of FIRST sources. A caveat in our bias 
analysis is that this bias may also depend on the actual 
angular extent of the unresolved sources. Therefore, the 
bias expectation value for a given flux may depends on 
a "hidden" parameter which value is not measured, and 
cannot be entirely removed from the data. 

10 Dominated by resolved sources. 

11 The fluxes in Tables [2] and [j] are not corrected for this bias. 
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ViRST^y] 

Fig. 4. — The FIRST peak specific flux vs. the same quantity 
from the NVSS catalog for the 4367 matched unresolved sources 
(see text). The solid line shows the 1:1 line, while the dashed lines 
represent 4tr below and above the 1:1 line. Open boxes shows the 
43 sources that are variable by more than 4<r, while all the other 
matched sources are marked by dots. We note that there are more 
variable sources below the lower line than above the upper line (see 
text) . 

The FIRST vs. NVSS peak flux-density measurements 
of the matched point sources are presented in Figure |3J 
Based on this plot we estimate the standard deviation, 
er, of the differences between FIRST and NVSS specific 
fluxes as a function of FIRST flux. This is done by cal- 
culating the 68-percentile range in the flux-flux plot as 
a function of the FIRST flux density, /first- We divide 
the 68 percentile by 2 to estimate the standard deviation, 
and then fit a first order polynomial to the logarithm of 
the standard deviation estimator as a function of flux 
density. We find that the relative errors associated with 
these difference measurements are well represented by 



a/f « 10 



-0. 6652-0. 0064/p 



(1) 



Here we define variables as objects for which the FIRST 
vs. NVSS flux difference is largeiFl than 4cr. We found 
43 such variable sources, which are listed in Table [3] We 
inspected the radio images of all these variable sources by 
eye, and comments on individual sources appear in this 
table. We note that the theoretical errors are smaller, 
and have different functional forms, than those implied 
by Equation [T] However, in order to avoid any possible 
uncertainties in the comparison between the two cata- 
logs, we used the empirical errors. 

Although we attempt to correct for the flux bias be- 
tween FIRST and NVSS sources, in Figure 0] there are 
more variable sources below the lower 4-cr line than above 
the upper 4-cr line. This systematic difference may be re- 
lated to the complexity of the bias between the FIRST 
and NVSS measurements, mentioned earlier. Effectively, 
this systematic bias induces errors in the number of "sig- 
mas" in which a source is variable. However, the ratio 
of number of sources above the upper 4-cr line to that 
below the lower 4-cr line, is consistent with an additional 

12 Assuming Gaussian noise, 4<r corresponds to probability of 
= 1/15,000 while the number of measurements in our experiment 
(number of epochs multiplied by the number of sources) is 8734. 
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Fig. 5. — The structure function of variable sources (black circles) 
and non-variable sources (gray circles) as a function of time lag. 
The measurements are displayed only for bins in which the number 
of flux differences is larger than four. The errors are based on the 
standard deviation in each bin. The horizontal dashed line marks 
the zero structure function. 

systematic shift in the flux ratio of « 0.35cr. Therefore, 
we conclude that since we used a relatively large vari- 
ability threshold of 4cr, most of our variable sources are 
probably real. 

4. THE STRUCTURE FUNCTION 

Next, we calculate the mean structure function for all 
the 43 sources which exceed the 4-cr variability thresh- 
old in Figure HI As a reference we also calculated the 
structure function for all the 3906 "non- variable" sources 
defined here as sources which variability is less than 2cr. 

The structure function was calculated according to the 
following scheme. For each pair of matched FIRST and 
NVSS measurements, we calculate the time difference 
between the FIRST (^ p IRST ) and NVSS (C p vss ) epochs: 
Here, i is the index of the pair. 



At,; = t 



FIRST 

cp 



NVSS 
cp 



We also calculate for each pair 



A/i//i 



/NVSS,i — fi 
fi ' 



(2) 



where /nvss.i and /firsts are the NVSS and FIRST 
specific fluxes of the i-th source, and fi = (/nvss,i + 
/first, i)/2. Next, the structure function and its error 
are estimated in bins of 500 days, between zero and 
3500 days, by calculating the mean and standard devi- 
ation of lS.fi/fi for all the pairs i in the appropriate bin. 

Figure [5] shows the structure function for the variable 
(black) and non variable (gray) sources. The structure 
function in this time range is consistent with being flat, 
with mean relative variability of about 35%. However, 
the value of 35% probably represents our sensitivity for 
variability rather than some physical variability level. 
The only physically interesting fact is the flatness of the 
structure function. 

5. DISCUSSION 

We present versions of the FIRST and NVSS catalogs 
that contain the mean epoch in which each source was ob- 
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TABLE 3 
Variable sources 



J2000.0 



USNO-B1 



2MASS 



ROSAT 



SDSS 



RA 


Dec 


/first 


/nvss 


At 


B 2 


R-2 


J 


H 


K Dist 


9 


r 






deg 


deg 


mjy 


mjy 


day 


mag 


mag 


mag 


mag 


mag " 


mag 


mag 




229.95525 


-5.90828 


19.40 ± 0.14 


10.0 ± 0.5 


1736.7 




















34.81875 


—4.64270 


5.72 ± 0.14 


3.0 ± 0.7 


1270.9 




















208.61233 


0.73462 


5.20 ± 0.18 


2.4 ± 0.5 


760.9 




















120.64970 


7.55780 


5.19 ± 0.15 


11.4 ± 1.5 


1758.9 




















193.96310 


8.74396 


7.90 ± 0.19 


3.7 ± 0.4 


1288.4 












24.70 


22.76 




0.30 


158.12037 


8.98541 


8.60 ± 0.15 


16.6 ± 0.6 


1719.3 


19.6 


18.5 








19.28 


19.28 


0.454 


0.45 


155.94935 


10.65521 


5.77 ± 0.15 


13.4 ± 1.6 


1698.6 












24.42 


22.33 






225.14947 


13.45289 


6.08 ± 0.14 


3.0 ± 0.4 


1271.5 




20.0 








23.01 


21.37 






114.80608 


18.03965 


7.26 ± 0.16 


3.3 ± 0.5 


1762.5 


17.1 


16.9 


16.14 


15.49 


14.86 


17.36 


17.12 




1.00 


201.32666 


19.46406 


7.04 ± 0.14 


12.7 ± 0.5 


1712.8 




















206.50009 


19.73055 


5.78 ± 0.15 


2.7 ± 0.4 


1702.0 












24.12 


23.18 






244.32375 


20.57935 


7.82 ± 0.14 


3.7 ± 0.4 


1320.1 












22.70 


22.29 






222.17540 


20.67543 


6.88 ± 0.15 


3.6 ± 0.5 


1322.1 




















247.57423 


23.33178 


9.41 ± 0.36 


2.8 ± 0.4 


286.9 


20.9 


19.7 








20.87 


20.52 




2.60 


212.99535 


23.89795 


9.53 ± 0.14 


5.0 ± 0.4 


284.4 












22.19 


21.96 






252.75117 


24.26111 


8.79 ± 0.14 


4.6 ± 0.4 


197.4 












22.86 


21.67 






190.23317 


24.33961 


5.80 ± 0.15 


2.5 ± 0.4 


284.2 




















182.98029 


25.36394 


5.48 ± 0.15 


2.8 ± 0.4 


213.8 


20.2 


20.5 








20.83 


20.66 




0.30 


247.48198 


27.72523 


7.28 ± 0.16 


3.8 ± 0.5 


183.4 




19.7 








22.40 


20.78 






238.73764 


29.95945 


5.08 ± 0.14 


9.2 ± 1.0 


—754.7 


18.9 


18.7 








18.89 


18.87 


0.855 


0.85 


203.08994 


29.99894 


16.62 ± 0.59 


8.9 ± 1.1 


—745.7 




















136.90323 


30.25872 


9.35 ± 0.15 


4.5 ± 0.5 


—286.1 












22.13 


21.99 






203.03348 


30.69115 


9.26 ± 1.35 


4.6 ± 0.5 


—735.7 




















169.75685 


36.87706 


9.33 ± 0.15 


18.5 ± 1.1 


—289.7 


16.3 


15.7 


15.98 


15.25 


14.60 


25.12 


19.25 






170.99205 


38.43646 


6.17 ± 0.15 


3.4 ± 0.4 


—260.7 






















OO.OU1 / 1 


l .ou it U. 14 


A A 4- n £ 
4.4 It U.O 


ZUZ.O 




















119.02080 


39.25484 


6.36 ±0.14 


3.0 ±0.4 


218.6 


19.2 


17.6 








20.45 


19.62 






262.32334 


44.77948 


5.98 ±0.14 


3.1 ±0.5 


700.7 




















229.90208 


44.97024 


5.12 ±0.15 


2.7 ±0.4 


710.8 












24.15 


23.91 




0.10 


134.99454 


45.87719 


30.74 ±0.15 


19.6 ±0.7 


1207.4 


19.2 


18.6 








18.86 


19.05 


0.440 


0.50 


220.45533 


50.71170 


5.06 ±0.16 


2.2 ±0.4 


753.7 




















262.65758 


52.24396 


7.32 ±0.14 


3.9 ±0.4 


766.6 






16.57 


15.86 


15.01 










211.65318 


52.68479 


8.42 ±0.22 


3.8 ±0.4 


772.9 




















159.71340 


53.49093 


8.53 ±0.15 


34.1 ± 1.8 


765.5 


19.9 








53.2 










226.56440 


53.92288 


5.30 ±0.15 


12.4 ±1.1 


770.9 


21.6 










21.75 


21.57 




2.40 


179.14175 


54.63838 


12.54 ±0.27 


7.2 ±0.5 


779.5 




















157.78099 


56.75153 


6.56 ±0.18 


3.3 ±0.4 


780.2 


20.6 


18.9 








20.84 


19.82 






157.96070 


56.85260 


6.29 ±0.20 


3.2 ±0.4 


780.2 




















160.55682 


59.35609 


5.51 ±0.15 


10.0 ±0.9 


2642.4 




















151.08074 


59.53517 


8.38 ±0.16 


3.5 ±0.4 


3075.2 




















185.45517 


61.59247 


7.51 ±0.15 


3.9 ±0.4 


2655.3 












23.10 


21.16 






248.86649 


62.64994 


14.23 ± 1.51 


7.1 ± 1.1 


2665.3 




















124.64543 


63.90501 


5.21 ±0.14 


2.2 ±0.4 


3192.3 





















Note. — List of 43 sources which vary by more than 4cr between the FIRST and NVSS epochs. The table is sorted by declination. Columns 
description: f is the peak specific flux and its error. The subscript indicate the catalog name. At is the time between the FIRST and NVSS 
observations. The position of each source was cross correlated with various catalogs, including the USNO-B1 (Monet al. 2003), 2MASS (Skrutskic 
ct al. 2006), ROSAT bright and faint source catalogs (Voges et al. 1999; Voges et al. 2000), and the Sloan Digital Sky Survey (York et al. 2000). 
In case counterparts are found we list their USNO-B1 E>2 and i?2 magnitudes, 2MASS J, H and K magnitudes, distance from ROSAT source, and 
SDSS g and r magnitudes and rcdshifts. We use search radius, relative to the FIRST catalog position, of 60 for ROSAT and 2.5 ' for all the other 
catalogs, z is the SDSS spectroscopic rcdshift of the source, while z p h is the photometric rcdshift of the source based on the SDSS colors. We use 
a photometric rcdshift estimator for quasars which is described in Ofek et al. (2002). The photometric rcdshift is calculated only if the source is 
indicated as a possible quasar in the SDSS database. Wc note that non of these sources are associated with a GB6 source (e.g., Gregory ct al. 1996) 
or a known pulsar. 
Comments on individual sources: 

RA— 155.94935 deg, Dec— +10.65521 deg : The peak flux wc measure in the NVSS image is a factor of two lower than the flux stated in the NVSS 
catalog, so this may be a constant source. 
RA= 247.57423 deg, Dec= +23.33178 deg : 
RA= 203.03348 deg, Dec= +30.69115 deg : 
RA= 159.71340 deg, Dec= +53.49093 deg : 
RA= 226.56440 deg, Dec= +53.92288 deg : 
RA= 157.78099 deg, Dec= +56.75153 deg : 
RA= 157.96070 deg, Dec= +56.85260 deg : 
RA- 124.64543 deg, Dec- +63.90501 deg : 



Near the strong source 3C 340 - NVSS and FIRST images are noisy. 
Radio images arc noisy. 

This is possibly a radio SN in the outskirts of NGC 3310 (Argo ct al. 2004). 
NVSS image shows a double source. 
Extended emission 15 arcmin from source. 
Extended emission 4 arcmin from source. 

FIRST flux may be influenced by sidelobes from 3C 343 (4.5 Jy, 9 arcmin to the NW). 
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served. We use these catalogs to look for variable sources, 
and we construct the structure function for these objects. 
We show that the structure function is flat on time scales 
between about half a year and five years. 

It is well known that the structure function of variable 
radio sources rises on time scales of days to tens of days 
(e.g., Qian et al. 1995; Gaensler & Hunstead 2000; Lovell 
et al. 2008; Ofek et al. 2011). Intrinsic variability of 
compact radio sources, which are mainly Active Galactic 
Nuclei (AGN), on days time scales would imply that the 
sources have small physical size. This in turn requires 
a very high rest-frame brightness temperature (TB lies t), 
orders of magnitude above ~ 10 12 K which is the limit 
for an incoherent synchrotron source (e.g., Kellermann & 
Pauliny-Toth 1969; Readhead 1994). Therefore, most or 
all of the variability of radio sources, below ~ 5 GHz, on 
these short time scales is presumably due to scintillations 
in the ISM. Moreover, based on causality arguments, in- 
trinsic variability of compact radio sources is expected 
only on time scales (ry) larger than 

( V ^ -3/2 ( 1 + Z N 3/2 ( T ^ -1/2 

x (io) (—) (ic^kJ day ' (3) 

where A/ is the variation amplitude in specific flux, c?i um 
is the luminosity distance (normalized at z = 1), T> is the 
Doppler factor of a relativistic motion in the source, and 
z is the source redshift. In Equation Truest is given 
in the rest frame and all the other parameters are in the 
observer reference frame. We note that Eq. [3] is derived 
from the relation for the brightness temperature, and 
replacing the source size by cry, where c is the speed of 
light. 

In contrast, observations of quasars and BL Lac ob- 
jects performed in the 4.8-14.5 GHz range showed that 
the structure function saturates only on time scales be- 
tween a year and ten years (Hughes, Aller & Aller 1992). 
Therefore, the fact that we do not see any significant 
rise in the structure function on time scales of months to 
years suggests that, at 1.4 GHz, the amplitude of intrin- 
sic variability relative to scintillations is small. Alterna- 
tively, it may suggest that the 1.4 GHz power spectrum 
of AGN radio variability is consistent with a white-noise 
power spectrum, rather than the red-noise power spec- 
trum typically seen at shorter wavelengths (e.g., Giveon 
et al. 1999; Markowitz et al. 2003). We note that 
Padrielli et al. (1987) reported on a class of sources (de- 
noted "C-BBV" in their terminology) which show cor- 
related intrinsic variability in low (0.4 GHz) and high 
(14.5 GHz) frequencies. However, these sources are a mi- 
nority among the variable sources in their sample. 

Figure |6] shows the fraction of variables as a function 
of Galactic latitude. The first bin contains two variables 
out of 42 objects (a fraction of 0.048). However, the ex- 
pectation value in this bin, estimated based on the mean 
fraction of variables (dashed line in Fig. [6] 0.0109) mul- 
tiplied by the number of sources in the first bin (42) is 
0.457. Assuming a binomial distribution, the cumulative 
probability to observe > 2 events, given an expectation 
value of 0.457, is 7.7%. This rules out the null hypothe- 
sis that the low-latitude variable-fraction is drawn from 
a uniform all-sky distribution at the 92.3% confidence. 
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Fig. 6. — The fraction of variables, relative to all the non- variable 
and variable sources, as a function of Galactic latitude. The er- 
ror bars represent the 1-cr errors on the fractions (Gehrels 1986). 
The dashed horizontal line marks the mean "all-sky" fraction of 
variables in our sample (i.e., 43/3949 = 0.0109). 



Therefore, a larger sample is required in order to confirm 
the earlier claims that the fraction of variables is larger at 
low Galactic latitude (e.g., Gaensler & Hunstead 2000). 
If this excess is real, then a plausible explanation is that 
it is due to ISM scintillations which are more prominent 
at low Galactic latitudes. However, we cannot rule out 
that some of this excess in variability is due to a popula- 
tion of Galactic variable sources as suggested by Becker 
et al. (2010). 

Finally, we use our dataset to look for correlation of 
the variability amplitude with redshift. There are sev- 
eral factors that can contribute to such a correlation. 
For example, scintillations depends on the source angu- 
lar size (larger amplitude for smaller sources), intrinsic 
source size evolution, broadening due to scattering, and 
maybe even scintillation in the intra-galactic medium 
(which may depends on the He re- ionization). Figure [7] 
shows Afi / fi as a function of redshift for all nine variable 
sources for which we have a redshift estimate (Table [3]) . 
The Spearman rank (Pearson) correlation coefficient be- 
tween the redshift and the relative variability amplitude 
is 0.59 (0.79). In order to estimate the significance of this 
correlation we conduct 10 4 bootstrap simulations (Efron 
1982; Efron & Tibshirani 1993). In each simulation we 
select, for each source redshift, a random variability am- 
plitude from the list of nine amplitudes. We find that 
the probability to get a Spearman rank correlation coef- 
ficient > 0.59 is about 5%. Therefore, unlike Lovell et 
al. (2008) we do not find any strong evidence for corre- 
lation between scintillations and redshift. However, our 
sample is considerably smaller than the one presented by 
Lovell et al. (2008). We note that eight out of the nine 
sources in Figure [7] are found above Galactic latitude of 
40 deg, and one source is at Galactic latitude of about 
18 deg. Removing the single low galactic latitude source 
(the source at z w 1 in Figure [7]) reduces the Spearman 
rank correlation to 0.42 and therefore does not change 
this results significantly. 
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Fig. 7. — The variability relative amplitude, A/i//;, as a function 
of redshift for all nine variable sources that are quasars or quasar 
candidates (Table|3jl. The upper x-axis shows the angular diameter 
distance corresponding to the redshift, assuming WMAP fifth year 
cosmological parameters (Komatsu et al. 2009). 
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